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Abstract 

In a recent paper we presented analytic expressions for the axis potential, the 
disk metric, and the surface mass density of the global solution to Einstein's 
field equations describing a rigidly rotating disk of dust. Here we add the 
complete solution in terms of ultraelliptic functions and quadratures. 
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a. Introduction. An infinitesimally thin disk of dust rotating uniformly around its sym- 
metry axis is the simplest model of a rotating self-gravitating body. Within Newton's theory 
of gravitation this model is described by the Maclaurin solution of the Laplace equation. In 
the present letter we show that the model of a uniformly rotating disk of dust allows for an 
explicit analytic solution in Einstein's theory of gravitation as well. This seems to be the 
first rigorous global solution of the rotating body problem (for perfect fluids with p = 0) in 
general relativity. 

An approximate solution of the problem was given by Bardeen and Wagoner [1], [2]. 
By applying the Inverse (Scattering) Method to the related boundary value problem of the 
Ernst equation we found the solution in terms of two linear integral equations, a 'small' 
and a 'big' one [3]. The analytic solution to the small integral equation in terms of the 
Weierstrass function led to explicit expressions for the Ernst potential on the symmetry 
axis, the disk metric, and the surface mass density [4]. In the following we will show that 
the complete solution of the problem may be represented, up to quadratures, in terms of 
ultraelliptic functions. 

b. The metric. The line element reads in Weyl-Lewis-Papapetrou-coordinates: 

ds 2 = e - 2U [e 2k (dp 2 + d( 2 ) + p 2 d V 2 } - e 2U (dt + ad<p) 2 . (1) 

(We use units where Newton's gravitational constant G as well as the velocity of light c are 
equal to 1.) The metric functions e 2U , e 2k , and a/p depend uniquely on the normalized 
coordinates p/po, C/Po, an d the parameter 

p = 2Q 2 p 2 e^, (2) 

where Q, po and Vo are the angular velocity, the coordinate radius, and the 'surface potential' 
V (p) = U(p — 0, £ = 0, p), respectively. The disk is characterized by ( = 0, < p < po- 

c. The Ernst potential. The (complex) Ernst potential is defined by 

f = e 2U + ib, (3) 

where the imaginary part b is related to the metric function a according to 

a P = pe- 4C/ 6 f ; a c = -pe~ 4U b >p (4) 

i. e., a may be calculated from a given Ernst potential (and its derivatives) by a quadrature. 
The same holds for the metric function e 2k . Hence, an axisymmetric stationary solution 
to the vacuum Einstein equations is sufficiently characterized by its Ernst potential. Note 
that our global solution is completely described by the vacuum solution outside the disk 
including the boundary data on the disk. 

The Ernst potential for the dust disk solution is given by the expression 



/ = exp { p 

with 
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W = W 1 W 2 , W 1 = y/{X - C/Po) 2 + (p/po) 2 , W 2 = ^/l+p 2 (l + X 2 ) 2 , (6) 
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h = ln(^/l + ^(l + X^ + Ml + ^ 2 )) 
Wl + /z 2 (l + X 2 ) 2 



/i /i 

The upper integration limits X a and Xb in the first two integrals are hyperelliptic functions 
of two variables u and v (i. e. 'ultraelliptic' functions) defined by 

fdX x rdx fxdx X [XdX 

X\ X2 X\ X2 

where the variables u and v are given by 

' hdX } hXdX 



u 



r hdX r hXdX 

= IwT- " = < 10 » 



The third integral in (5) as well as the integrals in (10) have to be taken along the imaginary 
axis in the complex X-plane with W\ and h fixed according to WNi < (for p, ( outside 
the disk) and = 0. 

Eqs. (9) have exactly the form of Jacobi's famous inversion problem [5], which was solved 
by Gopel [6] and Rosenhain [7] who were able to express X a (u,v) and Xj,(u,v) in terms of 
(ultraelliptic) theta functions. In our case the moduli of the theta functions depend on p/po, 
C/po, and p. 

It should be noted that the integrations from Xi to X a resp. X 2 to X^ in (5) and (9) have 
to be performed along the same curve in the two-sheeted Riemann surface associated with 
W(X). This leads to a unique Ernst potential f(p/po, C/Po, A*)- More technical details and 
final expressions for the metric in terms of theta functions will be presented in a subsequent 
paper. 

d. The Newtonian limit /i < 1. At p, — 0, the solution may be expanded into a power 
series in p}l 2 : 

00 

/ = 1 + £ iWi (n+1)/2 . (11) 

n=l 

The coefficients f n {p/ Po,(/ Po) turn out to be elementary functions. A particularly simple 
result is obtained for fi to / 4 . From Eqs. (5) to (10) we get 

X a -X 1 = 0{ i ?l% X b -X 2 = 0(p^ 2 ), (12) 

and 

f^ e J-,) ^- X ^-^ dX + OiA. (13) 
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Expanding h, Xi and X 2 one can easily determine fi to / 4 . (Note that X = 0(1) in this 
integral.) The result for fi is 



A -l/fi±g"--I{i«,r. €+ [ e -( f + I,^ e ](i-, 1 rt} (u, 



—1 

with elliptic coordinates £, 77 defined by 



P = Po V /l + £ 2 V 1 - r 7 2 ' C = Po^, 0<^<oo, -1 < 77 < 1. (15) 
Now, the Newtonian limit Umi can be read off from the first order relation 

9u = -*/ = -(1 + nh) = -(1 + ^). (16) 

cr 

Obviously, Umi = c 2 /z/i/2 is the Maclaurin solution, i. e. the Newtonian solution of the 
problem of a rigidly rotating disk of dust. Note that in this limit 

2U M i(p = 0,C = 0) _ ^ 2 pI ( . 

»- C 2 - C 2 ■ y u ) 

(For the sake of clarity we have reintroduced the velocity of light c in the last two equations.) 

e. Discussion. In the parameter range < /i < /io = 4.62966184 ... [4] the solution 
is asymptotically flat and regular everywhere outside the disk. Its baryonic surface mass 
density is positive and vanishes at the rim of the disk. For fi fi (il finite) and p 2 + ( 2 ^ 
it approaches exactly the extreme Kerr solution (mass M — 1/ 2Q). The meaning of the 
solution for fi > fi deserves further discussion. 

The analytic expressions above allow for a detailed investigation of characteristic prop- 
erties, as for instance the formation of toroidal ergospheres, dragging effects, motion of test 
particles, and stability. 
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